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Q I Abstract 

Let A and i? be positive semidefinite matrices. We investigate the conditions under 
which the Lieb-Thirring inequahty can be extended to singular values. That is, for 
. which values of p does the majorisation a{B^AP) -<yj a{{BA)''') hold, and for which 

I ! values its reversed inequality a{BPAP) >-y; a{{BA)P). 
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■ The famous Lieb-Thirring inequahty |6j states that for positive semidefinite 
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matrices A and B, and p > I, Tt{AB)p < Tr{APBP), while for < p < 1 
the inequahty is reversed. Many generahsations of this inequahty exist [2117] . 
one of the most noteable being the Araki-Lieb-Thirring inequahty [1]. For 
positive matrices A and B, and any unitarily invariant norm ||| ■ |||, the fol- 
lowing holds (see also Theorem IX.2.10 in [3j): \\\{BAB)p\\\ < \\\BpApBp\\\ 
when p > 1, and the reversed inequality when < p < 1. This inequality 
can be equivalently expressed as the weak majorisation relation between sin- 
^ : gular values a{{BAB)P) a{BPAPBP). Here, a{X) a{Y) if and only 

5^ ; if J2^j=iaj{X) < Y.j=i(^j{y), ioT 1 < k < d, where aj{X) denotes the j-th 

largest singular value of X. 

In this paper we study the related question whether a majorisation relation 
exists between the singular values of the non-symmetric products BpAp and 
[BA)P. The latter expression is well-defined because the eigenvalues of a prod- 
uct of positive semidefinite matrices are real and non-negative. Our main result 
is the following: 

Theorem 1 Let A, B > be d x d matrices. For < p < 1/2, 

a{BPAn aiiBAY). (1) 
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In addition, if d = 2, the range of validity extends to < p < 1. 
For p > d — 1 and for p G Nq, the reversed inequality holds: 

a{BPAn aUBAy). (2) 



In the first half of the paper, we prove this Theorem for p satisfying the 
condition 1/p & Nq or 1/p > d — 1 or p & Nq or p > d — 1. We do so by 
chaining together two majorisations; in terms of the first inequality, ([T]), we 
chain together a{APBP) aP{AB) and aP{AB) a{{AB)P). While the first 
majorisation indeed holds generally and is a straightforward consequence of 
the original Lieb-Thirring inequality, see Theorem [2l the second majorisation 
turns out to be subject to the rather surprising condition on p (Theorem [3]). 
In the second half of this paper, we follow a different route and obtain validity 
of (P for < p < 1/2. 

Henceforth, we abbreviate the term positive semi definite as PSD. 
The following Theorem is already well-known: 
Theorem 2 For A, B PSD, andO <p <1, 

a{APBP) a^iAB). 
For p > 1, the direction of the majorisation is reversed. 

Proof. We only have to prove the statement for cti, i.e. the infinity norm ||.||. 
From that we can derive the full majorisation statement by using the well- 
known trick, due to Weyl, of replacing X by its antisymmetric tensor powers, 
as in yy. 

Consider < p < 1. By the Araki-Lieb-Thirring inequality for the infinity 
norm || ■ ||, we have 

\\AB^A\\P = \\{AB^AY\\ > \\APB^PAP\\. 

Noting that ||XX*|| = this gives, indeed, 

\\AB\\P > \\APBP\\. 

This inequality was first proven by Heinz (see Theorem IX.2.3 in [3]). For 
p > 1, the direction of the inequalities is reversed. □ 

For the second majorisation we need a lemma, which relates the question to 
a result by FitzGerald and Horn. 
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Lemma 1 Let {Xi)i be a sequence ofd non-negative numbers. The dxd matrix 
C with entries 

1 — AjAj 

IS PSD z/a G No or a > rf - 1. 

Proof. This expression can be represented in integral form as [5] 

1 

Ci,j = a J dt{t+{l-t)XiXj)''-\ 



Thus C is PSD if the integrand is. Since for < t < 1 the matrix (t + (1 — 
t)XiXj)ij is PSD and has non-negative entries, C being PSD follows from a 
Theorem of FitzGerald and Horn [5J that states that the g-th entrywise power 
of an entrywise non-negative PSD matrix is again PSD, provided either g G 
or q > d — 2. Here, q = a — 1, hence the condition is a G Nq or a > c? — 1. □ 

Theorem 3 Let X be a d x d matrix with non-negative real eigenvalues. For 
p in the range < p < 1, the majorisation 

holds, provided 1/p G Nq or 1/p > d — 1. 

For the range p > 1, the direction of the majorisation is reversed, and the 
conditions for validity are ^gNq or p > d — 1. 

Proof. Consider the case < p < 1 first. 

Again, we consider the inequality cri{X) < ai{XP), from which the majorisa- 
tion of the Theorem follows by the Weyl trick. 

An equivalent statement of the inequality is: \\Xp\\ = 1 implies \\X\\ < 1 
(obtainable via rescaling X). 

If we impose that X be diagonalisable, it has an eigenvalue decomposition 
X = SAS^^, where S is invertible, and A is diagonal, with diagonal entries 
Xk > 0. Then 



\\xp\ \ = i^{xpy{xp) < 11 

^APS*SA'' < S*S. 

Let us introduce the matrix A = S*S, which of course is positive definite, by 
invertibility of S. Thus the statement | \Xp\ \ = 1 is equivalent with A^AA^ < A. 
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Likewise, the statement ||X|| = 1 is equivalent with AAA < A. We therefore 
have to prove the imphcation 

A^AA^ <A^ AAA < A. (3) 



Now note that, since A is diagonal, the condition A^AA^ < A can be written 

as 

A' :=Ao(l-A?A^)f,=i>0, 
where o denotes the Hadamard product. Likewise, AAA < A can be written 
as 

A o (1 - A.A,)t=i > 0. 

In terms of A', this reads 

A'oC>0, 

with 



C :-- 



XiXj \ 



1 — Af A^ / • • 1 



Thus, by Schur's Theorem |4j, the implication ([3]) would follow from non- 
negativity of the matrix C. Using Lemma [1], we find that a sufficient condition 
is 1/p G No or 1/p > c/ - 1. 

Using a standard continuity argument, we can now remove the restriction that 
X be diagonalisable. 

The case p > 1 is treated in a completely similar way, but relying instead on 
the non-negativity of the matrix 

/1-AfA^y 



XiXj 



□ 



For all other values of p than the mentioned ones, the matrix C encountered in 
the proof is in general no longer non-negative. Likewise, for these other values 
of p, counterexamples can be found to the inequality that we wanted to prove 
here, so the given conditions on p are the best possible. 

Combining Theorem 2 and Theorem 3 immediately proves Theorem 1 for 
1/p G No or 1/p > d — 1 or p E No or p > d — 1. 



To prove the remaining case covered by Theorem [H we derive several equiva- 
lent forms of the inequalities ([1]) and ([2]). We again only need to treat the ai 
case, as the full statement follows from it using the Weyl trick. 
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Consider first the case < p < 1. Then we need to consider 

\\BPAP\\<\\{BAY\\, (4) 

since the largest singular value is just the operator norm. 

As a first step, we reduce the expressions in such a way that only positive 
matrices appear with a fractional power. 

By exploiting the relation ||X|| = (|4j) is equivalent to 

\\APB^PAf\ \ < \ \{ABy{BAf\\, 

which, by homogeneity of both sides, can be reformulated as 

\\{ABY{BAY\ \ < 1 =^ WAPB^PA^W < 1, 

and, in terms of the PSD ordering, 

{ABf{BAY < 11 ^ A^B^^AP < 11. (5) 

Lemma 2 For any A > and B > 0, there exist diagonal A > and in- 
vertible S such that A = SS* and AB = SAS^^, and, consequently, B = 

S-*AS-\ 

Proof. Let AB = TAT^^ be an eigenvalue decomposition of AB. Because 
A and B are PSD, the eigenvalues of AB are non-negative, hence A > 0. 
Assuming that all eigenvalues of AB are distinct, we show that T^'^AT^* is 
necessarily diagonal. 

Indeed, from AB = TAT'^ follows T'^AT-* T*BT = A. The factors X = 
T~^y4T~* and Y = T*BT are positive definite, and positive semidefinite, 
respectively, since they are related to A and i? by a ^-conjugation. Now note 
that A is diagonal and all its diagonal elements are distinct. This implies 
that X and Y, both Hermitian, are themselves diagonal. This follows from 
taking the hermitian conjugate of XY = A, YX = A, and noting that the 
two equations taken together imply that X and Y commute and are therefore 
diagonalised by the same unitary conjugation. Then we see that the product 
XY must also be diagonalised by that same unitary conjugation. However, 
XY = A is already diagonal, so that X and Y must be diagonal too. 

By a continuity argument, we see that there must exist a T diagonalising both 
AB (via a similarity) and A (via a *-conjugation) even when the eigenvalues 
of AB are not distinct. 

The lemma now follows by putting S = TX^/^. □ 
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Using the Lemma, the left-hand side (Ihs) of ([5]) can be rewritten as 



{ABy{BAy = {SAs-y{s-*As*y 

= SAPS-^S-*APS*. 
The condition {ABy{BAy < 11 then becomes 

APS-^S-*AP < S-^S-*, 

which turns into 

A^CA^ < C 

on defining C = S'^S'* > 0. 

Similarly, the condition of the right-hand side (rhs) of A^B^^A^ < H, can 
be rewritten as B^^ < A~^p, or 

{S-*AS-^yp < {SS*)-^P = {S-*S-^yp. (6) 

Using the polar decomposition, we can put 5""* = f/C^/^, where f/ is a uni- 
tary matrix. Then the condition of the rhs becomes ([/C^/^AC^/^f/*)^^ < 
{UCU*yP, or 

{C^/^AC^/yP < C^P. (7) 

Thus, implication is equivalent to 

A^CA*' < C ^ {C^/^AC^/yP < C^P, (8) 

for < p < 1, and C > 0, A > 0. 

On left- and right- multiplying both sides of the Ihs of ([8]) with C^/^, we get 

{C'/^A^C'/y <C^^ (C'/^AC'/^P < C^P. 

By putting A = C^l"^ and B = A^, this becomes 

[ABAf <A^^ {AB^/PAyp < A^p. 

In this equivalent form, it is now easy to prove ([T]) for p < 1/2. 

Proof of Theorem 1 for < p < 1/2: By operator monotonicity of the square 
root, {ABAy < A^ implies ABA < A^. Dividing out A on both sides, this 
is equivalent with B < 11. This implies B^^^ < 11, for all p > 0, and thus 
AB^/PA < A^. S ince < p < 1/2, operator monotonicity of the 2p-th power 
finally implies {AB^'PAfP < A'^p. □ 
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For d > 2 and l/2<j9<l,we have found counterexamples. To narrow down 
the search for counterexamples, we semi-intelligently chose a random positive 
diagonal d x d matrix D and a random (i- dimensional vector ip to construct 
A and B matrices: 



The condition {ABAf < A'^ is equivalent with \ \A-^BA^BA-^\\ < 1 and is 
thus satisfied by construction. However, with high probability A and B are 
found that violate [AB^/^ A^'p < A^'p. As the violations are extremely small, 
all calculations have to be done in high-precision arithmetic (we used 60 digits 
of precision) 0- This numerical procedure yielded counterexamples for d = ?> 
and p between 0.89 and 1. 

In a similar way counterexamples can be found in the regime d > 2 and p > 1. 
For p > 1 , we find by a similar reasoning that the reversed inequality of (jlj) is 
equivalent to the converse of ([5]), and therefore to the converse implication 



For d = 3 we have found counterexamples up to p = 1.25, but no higher. It 
is therefore imaginable that the second majorisation inequality in Theorem 1 
could be valid under more general conditions, e.g. for p > 2 perhaps. For the 
time being, this problem is still open. 
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